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We consider the braiding of Kramers pairs of Majorana bound states. We derive the most gen¬ 
eral transformation on the many-body ground state that is applied as the result of such a braiding 
process. The result is derived in the context of a simple toy model, but we will show that it has 
the most general form that is compatible with local and global conservation of electron parity. In 
accordance with earlier work the resulting transformation turns out to be path dependent, which 
shows that Kramers pairs of Majorana bound states cannot be used for topological quantum com¬ 
putation. We also discuss under which conditions the result is path independent and corresponds 
to two independent exchanges of pairs of Majorana bound states. 


For certain classes of topological superconductors (D, 
BDI), Majorana bound states (MBS) appear at topolog¬ 
ical phase boundaries or in vortices. Those MBS are a 
signature of the non-trivial topological phases and they 
can be used for topological quantum computation |Tj. As 
such MBS have received a lot of attention, both from the 
theoretical and the experimental side nu. 

Another closely related topological class is DHL This 
class has time reversal symmetry that squares to minus 
one and exhibits Kramers pairs of MBS at topological 
phase boundaries or in vortices. Even though there have 
not been experimental attempts yet to realize this topo¬ 
logical phase, there have been many theoretical proposals 
on how to obtain this phase [K tiTfil . A natural question to 
ask is whether those Kramers pairs of MBS can be used 
for topological quantum computation similar to MBS ap¬ 
pearing in classes D and BDI. This question is equiv¬ 
alent to the questions whether Kramers pairs of MBS 
are non-Abelian anyons and whether their adiabatic ex¬ 
change generates a transformation that is independent of 
the details of the exchange path. 

In earlier work we argued that this cannot be the case 
because Kramers pairs of MBS have a local degree of 
freedom that can be manipulated adiabaticallyp3|. This 
makes any exchange of two such Kramers pairs path de¬ 
pendent, which means that they are not anyons or equiv¬ 
alently that they cannot be used for topological quantum 
computation. 

In this article, we study the general transformation 
that occurs when braiding two Kramers pairs of MBS. We 
study it in the context of a simple toy model, but this will 
not limit the generality of the results because even in this 
simple model we find the most general transformation 
that is allowed by local parity conservation. 

The toy model we study is the time-reversal-symmetric 
analog of braiding by tuning couplings m • It consists 
out of four Kramers pairs of Majorana fermions \ii Xi = 
TxiT^ 1 , which we write as a vector Xi = {XiiXi) T ■ The 
most general time-reversal-symmetric coupling between 
two Majorana Kramers pairs is 

Hj k = HxJ^ze^^Xk- (1) 



FIG. 1. Each circle denotes a site with a Kramers pair of 
Majorana fermions. Lines between the sites denote couplings. 
As a single coupling is switched from position to the next, the 
corresponding zero energy Kramers pair (color) moves the 
opposite direction. Note that the sites label the static basis 
and therefore the numbers don’t change. 


In the toy model there are only three coupling terms, H, n, 
F ?42 and i? 43 - In order to have the appropriate ground 
state degeneracy at all times, there is at least one non¬ 
zero coupling and at most two at each point in parameter 
space. Furthermore, we take the Hamiltonian with two 
couplings to be of the form 

H = cos OHij + sin OHik, (2) 

such that 0 controls the relative strength of the couplings. 
The effect of switching from one coupling to another, by 
means of tuning 0 from 0 to 7r/2, is to move a zero energy 
Kramers pair from one site to another. The situation is 
depicted in Fig. [T] 

In order to study the adiabatic braiding process we 
have to specify a basis at each point in parameter space. 
At each such point we have two zero energy Kramers 
pairs of Majorana bound states, which we denote with 
X\ and X' 2 - We assume that those are local and initially 
start out as.X’i~xi,-^ 2 ~X 2 i where ~ means up to 
(generally different) local rotations e laiJy . We construct a 
local basis by forming local fermions D v = |( X^ +iX v ), 
where we used a Greek index which is either 1 or 2 to 
distinguish it from the Latin indices which run from 1 to 
4. The basis we use at each point in parameter space is 


2 


now given as 


|00> 

| 01 }=£> 1 | 00 ), 

H)) = I)\ 00) : 

|ll)=^Dt|00). 

Note that there is a problem with our local basis choice 
([3]) when completing one braid. By definition our X q 
are local and after the braiding is done they’ll end up 
as Xi ~ X 2 and X 2 ~ Xi- This means that our local 
basis ([3]) does not go back to itself after completing a 
loop in parameter space. We can correct for that by 
an explicit basis transformation at the end. This way 
the unitary transformation, U, due to braiding, factors 
into local Berry phases, Ui oca .i, followed by a global basis 
transformation, B. We have 

U = BU local . (4) 

In principle one could perform a basis transformation 
that exchanges X\ and X 2 at an earlier point during 
the exchange. In that case one still has to make sure 
that the phases of the states at the end match up exactly 
with the ones in the beginning, otherwise one needs to 
do another basis transformation. Therefore it is most 
convenient to perform everything as a single bases trans¬ 
formation, which includes potential phases, at the end of 
the process. 

In our earlier work m we studied Ui oca \. It generally 
takes the form 

U local = e^ x ^e^ x ^, (5) 

where the piS are calculated as 

Vv = \j{X r] ,\7 x X v }- dA, ( 6 ) 
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FIG. 2. Graphical representation of Q. A loop in parameter 
space is achieved by a continuous exchange followed by a basis 
transformation. The circles represent Kramers pairs of Ma- 
jorana fermions. The ambiguity of dividing a circle into two 
represents the gauge freedom in choosing Kramers partners. 
Two different lines and colors are used to divide the circles to 
distinguish the two Kramers pairs. 


A pictorial representation of the decomposition Q is 
given in figure [2j where the circles represent Kramers 
pairs and the gauge freedom in the choice of Kramers 
partners is represented by choice of angle at which the 
circles are split into two halves. 

One might be worried about the transformation prop¬ 
erties of equation 0 under gauge transformations. First 
of all it is clear that the decomposition can only hold as 
long as we restrict ourselves to local gauge transforma¬ 
tions. Under a gauge transformation W( A) the operator 
U should transform as 

U ^W{X initial )UW\\ 

initial ). (8) 


with C being the path in parameter space during the 
exchange. We refer to this as mixing by an angle tp 1 or 
tp 2 of the Kramers pair 1 and 2 respectively. Note that 
the integration path is not closed. Therefore a gauge 
transformation will change ([6]) according to the initial 
and final gauge choice. We will discuss this issue in detail 
later. In this work, we will focus on finding the basis 
transformation B in order to calculate Q. 

To find the basis transformation we have to determine 
the final form of X q . The transformation B is then de¬ 
termined as the transformation that fulfills 

= w; nitial . (7) 

This way B is, of course, only defined up to an abelian 
phase, but this phase will be strongly system dependent 
in any case and we do not consider it. From the definition 
it is also obvious that B depends on the gauge choice 
at the initial and final position. We will discuss this 
dependence in detail below. 


which means that the transformation U is basis depen¬ 
dent as one would expect. We explicitly show that this re¬ 
lation holds for our decomposition ©j under local gauge 
transformations, which have the form 

W local = eBM x ^ e BW x ^. (9) 

This is important, because we will show below that B and 
Uiocai do not satisfy the same transformation property 
individually. 

We denote Wi ocal (A initial/final ) by W i/f respectively. 
It follows directly from BX^ nal B^ = X“ litial 
that Wj B W}W{ X® nal w/W f B^W} = W ; X™ tial W i t . 
Hence B transforms as B —► WiBWj. If one replaces 

f Y„ X„ by W local X v Wl cal , Wl cal X, W loca i in © one eas¬ 
ily finds that the local phases transform as + 

(2/r,(A fina i) - 2/ r) (A ini ti a i)), from which it follows that 
Uiocai transforms as C/ loca i —► W f U\ oca \W^. Hence nei¬ 
ther B nor C/iocai transform according to © , in particular 
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because they depend on the gauge choice at Afj na i. This 
dependence however drops out of their product, such that 
U indeed fulfills ([8]), proving that our decomposition 0 
is valid for any local gauge choice. 

Before we continue and calculate B let us comment on 
a particular interesting local gauge choice. This gauge 
choice is defined by U\ oca \ = 1, which is always possible. 
This way one simply has U = B, which is useful for nu¬ 
merical calculations, because in this way the calculation 
of the transformation does not require taking derivatives. 
A discretization of a path in parameter space only needs 
to be fine grained enough in order to fix the gauge choice. 

To calculate B we have to find X)) nal and therefore we 
need to understand how X v changes during one switching 
process 0 when taking 9 = 0 —> . We require local¬ 
ity of the zero energy Kramers pairs. Therefore one of 
the X v will stay constant through the switching process 
(up to parameter dependent local gauge choices). The 
other of the X n will move between sites and by solving 
[H, X,,] = 0 we find 

X v = e iaa « (cos de^^Xj ~ sin 9e i0ik ^Xk) , (10) 


where a is an arbitrary gauge choice that may depend on 
the /3 s and 9. 

We want to patch three switching processes together 
such that X 7J is continuous. In order to do that conve¬ 
niently we pick a particular parameter dependence for a, 
such that X r/ (9 = 0) = ±Xj and X v (9 = |) = q=Xfc- 
The sign change is motivated by analogy to switching 
processes of single Majorana fermions in class D systems. 
A particular gauge choice of a that gives the desired X v 
at 9 = 0,, f is a± = —f3ij cos 9 — fin- sin 9 + (1 =F 1) § ■ 

Table ITT shows X v initially and after each switching 
process. It has exactly the same structure as for single, 
class D, MBS. As the transformation from initial to final 
Kramers pairs of MBS we get 


X initial _ _, _ 

i — Xi -t — X2, 

X initial _ _. v final _ 

2 — Xi ~> — ~Xl- 

Therefore the basis transformation that generates 
transformation is simply 


( 11 ) 

this 


B _ e f (X1X2+X1X2) 


( 12 ) 


which is structurally the same as two independent braid¬ 
ing transformations of pairs of MBS . According to 0 
the total braiding transformation is 


JJ _ gf (XlX2+XlX2) e ^XlXl + ¥x2X2_ 


(13) 


The question arises how general this transformation is, 
because we obtained it in the context of our simple toy 
model. From the construction it is clear that we can al¬ 
ways choose a parameter dependent basis such that Eqs. 
(111 and (12) are true. Therefore the question reduces to 
whether U\ oca \ can be more complicated than ([ 5 ]). Since 



#43 

#41 

#42 

#43 

Xi 

XI 

~X3 

~X3 

X2 

x 2 

X2 

X2 

-Xi 

-Xi 


TABLE I. Instantaneous X for a braiding process requiring 
them to be continuous. The process is of the form that is 
depicted in figure [l] 


we assume that both Kramers pairs of Majorana fermions 
are always decoupled, the local parity operators iX v X v 
are conserved quantities. The only local transformations 
which we can construct out of those are given by U\ oc& \. 
Hence our toy model already describes the most gen¬ 
eral possible braiding transformation for Kramers pairs 
of MBS. 

Because the phases in Eq. ([5j) are path dependent, 
the exchange of Kramers pairs of MBS cannot be used 
for topological quantum computation. The natural ques¬ 
tion arises whether there are additional conditions un¬ 
der which the phases become path independent. We will 
show that this is the case in the absence of local mixing. 
By local mixing we mean the transformations that can 
arise when changing the parameters along any closed tra¬ 
jectory in parameter space that does not braid Xi,X 2 . 
The transformation describing this process will simply be 
U = Hiocap with the difference that the integral in Eq. 
([6]) is over a closed trajectory in this case. The absence 
of local mixing means that is zero for any such closed 
loop. This result is useful because the presence or ab¬ 
sence of local mixing is a property of individual Kramers 
pairs of MBS, but we can use it to make a statement 
about braiding of several Kramers pairs. 

Before we discuss the implications of the presence or 
absence of local mixing, note that if the phases in Eq. ([ 5 ]) 
are path independent, then there cannot be local mixing, 
because otherwise we could always add a closed path with 
non-vanishing mixing phase to the open path in Eq. 0. 
violating the path independence of the phase. The ab¬ 
sence of local mixing is therefore necessary in order to 
have constant phases in Eq. We now proceed to 

show that it is also sufficient and moreover that in the 
absence of local mixing the braiding transformation can 
be reduced to two independent braidings of xi 1 X 2 and 
Xi,X2 respectively. 

We first show that are path independent. Assume 
we have two braiding paths 1 and 2 such that the mix¬ 
ing angles along those paths are ip v and (p v respectively. 
The braiding transformations are given by U and U. We 
can now form a local mixing operation U'U with mix¬ 
ing angles ip v — (p v , but since there is no local mixing by 
assumption we have ip v — (p v = 0. Therefore the are 
path independent. The situation is illustrated in figure 

©. 

We now argue that (p 1 = — ip 2 - In our simple toy model 
this can be checked by means of a straightforward calcu- 
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FIG. 3. (a) Two exchanges along different paths can be com¬ 
bined into two local operation (solid and dashed line). Conse¬ 
quently the mixing angles are path independent if local mixing 
angles are zero, (b) One Kramers pair gets moved out of the 
way for the other to make a loop. During the first half pair 1 
rotates by tp\ and during the second half it rotates by ip2 ■ If 
there is no local mixing this implies tpi + ip2 = 0. 


lation of t/i oca i, but we will give a more general argument 
instead. In the way we decomposed the braiding trans¬ 
formation, {/i oca i does not depend one the presence of the 
second Kramers pair of MBS. We can therefore imagine 
moving the second Kramers pair a little out of the way, 
such that the first Kramers pair can make a loop, which 
does not encircle the second one, but passes through the 
initial position of the second one. During the first half of 
this loop Xi will acquire a mixing angle which is the 
same as the mixing angle it acquires during an exchange. 
During the second half of this loop X\ will acquire a mix¬ 
ing angle y> 2 , which is the same as the mixing angle X 2 
acquires during an exchange. In the absence of local mix¬ 
ing the overall mixing angle for the loop is <£1 + <£2 = 0. 
This situation is illustrated in Fig. [3j 

We showed that in the absence of local mixing ipi = 
—<£2 in which case we can bring Eq. ( |13[ ) to the form 
of Eq. (121 by means of the local basis transformation 


W = e 4 XlXl * X 2 X 2 ^ suc h that 

WUW * = e ? (xix2+xix2)_ 


(14) 


This is only meaningful if the ip„ and therefore W do 
not depend on the specific braiding path, which we also 
showed. The question arises what the physical mean¬ 
ing of this particular basis choice is. This question can 
only be answered meaningfully if we also specify a con¬ 
dition that ensures the absence of local mixing. In our 
earlier work we specified sufficient symmetry condi¬ 
tions that guaranteed the absence of local mixing and in 
this case the basis choice leading to Eq. (141 is the one 


respecting this additional symmetry. 

It should be noted that a symmetry condition (or some 
other kind of condition) that guarantees the absence of lo¬ 
cal mixing does not need to be global. It only needs to be 
true locally around the Kramers pair. That means that 
one could try to engineer a system such that the sym¬ 
metry condition is approximately satisfied close to the 
topological phase boundaries, where the Kramers pairs 
of MBS are localized. If this remains true as the phase 
boundary and hence the Kramers pair moves, than one 
can get path independent braiding statistics. 

Even with the ideal transformation of Eq. (14) the 
question remains whether the braiding can yield some¬ 
thing more interesting then moving localized fermions 
around. The answer is yes, but it would require schemes 
that initialize and read out non-local states. One can for 
instance think of bringing two Kramers pairs of MBSs 
close together such that they split in energy and one 
can initialize them in the lowest energy state of the split 
states, after which one brings them apart again. Oper¬ 
ations on those states can then be performed similarly 
to the operations performed on class D Majorana qubit 
systems. The non-local states in class Dill systems are 
actually entangled with respect to the local bases. How¬ 
ever, local mixing will decohere the superpositions and 
therefore absence of local mixing is crucial for any oper¬ 
ations involving entanglement in the local basis. 

In conclusion, we derived the most general braiding 
transformation for Kramers pairs of Majorana bound 
states. This transformation naturally decomposes into a 
transformation that describes the independent exchange 
of two pairs of MBS and additional local rotations of the 
Kramers pairs that are being exchanged. The angles of 
those local rotations are generally path dependent and 
independent for the two Kramers pairs. The fact that 
the angles are path dependent shows again that Kramers 
pairs of MBS cannot be used for topological quantum 
computation. The derived form of the transformation is 
the most general one allowed by global and local par¬ 
ity conservation and therefore the result is not limited 
to the simple toy model, which we used to study it. We 
also defined the concept of local mixing, and showed that 
whether or not the phases are path dependent is equiva¬ 
lent to the presence or absence of local mixing. We also 
argued that in the absence of local mixing the exchange 
transformation reduces to two independent exchange of 
pairs of MBS without any path dependence. This is help¬ 
ful because it might be easier to check/engineer the pres¬ 
ence or absence of local mixing, because it is a local prop¬ 
erty. 
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